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Abstract 

The migration of Siamese mud carp (Henicorhynchus siamensis and H. lobatus), two of the most economically important fish 
species in the Mekong River, was studied using an otolith microchemistry technique. Fish and river water samples were 
collected in seven regions throughout the whole basin in Thailand, Laos and Cambodia over a 4 year study period. There 
was coherence between the elements in the ambient water and on the surface of the otoliths, with strontium (Sr) and 
barium (Ba) showing the strongest correlation. The partition coefficients were 0.409-0.496 for Sr and 0.055 for Ba. Otolith Sr- 
Ba profiles indicated extensive synchronized migrations with similar natal origins among individuals within the same region. 
H. siamensis movement has been severely suppressed in a tributary system where a series of irrigation dams has blocked 
their migration. H. lobatus collected both below and above the Khone Falls in the mainstream IVlekong exhibited statistically 
different otolith surface elemental signatures but similar core elemental signatures. This result suggests a population 
originating from a single natal origin but bypassing the waterfalls through a passable side channel where a major 
hydroelectric dam is planned. The potential effects of damming in the IVlekong River are discussed. 
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Introduction 

Migratory behaviors are extremely common in fisheries, and 
understanding fish migrations is vital to sustainable management 
and conservation of fisheries resources. While diadromous species 
and populations have long been studied mainly in temperate 
regions, potamodromous migrations that are common in large 
tropical rivers, such as the Amazon and Mekong Rivers are much 
less understood [1-4]. 

The Mekong River, an international river rising in China and 
flowing 4,200 km through five other countries of Southeast Asia, is 
renowned as the world's largest inland fishery, as well as for 
harboring the second most diverse fish fauna after the Amazon 
River [5] . It produces over 2 million tonnes of wild fish and other 
aquatic animals each year within its catchment [6] and supports 
more than 1,200 fish species, although many have not yet been 
described [7]. 

There is a clear and urgent need to better understand fish 
movements for two main reasons. Firstly, a significant proportion 
of Mekong fishes are migratory, specifically potamodromous (i.e.. 



migratory wholly in freshwater) [8], and, more importantly, this 
enormous fisheries production relies heavily on the migratory 
habit of these species [9], [10]. They are harvested mostly with 
gears and techniques developed to capture fishes during then- 
seasonal migrations [11-14]. Secondly, the Mekong fishes and 
fisheries now face an imminent threat from massive hydropower 
development. Currently blocked by mainstream dams only in the 
uppermost Chinese territory, the Mekong River has 1 1 more dam 
proposals on the mainstream and more than 100 in tributary 
systems in the lower basin, all of which are scheduled to be 
installed within the next decade [15]. These dams will impact 
fishes directiy as a barrier to their migrations or indirectiy through 
creation of huge impoundments and alterations to the natural flow 
regime, leading to a substantial loss of fish species and their 
habitats [10], [16-18]. 

Of the vast multitude of fish species in the Mekong River, two 
closely-related species, Henicorhynchus siamensis and H. lobatus, 
are of special concern given the rapid rate of hydropower 
development. Inhabiting the Mekong and Chao Phraya basins [H. 
siamensis) and the Mekong basin {H. lobatus), these two small- 
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Figure 1. Fish sampling sites in the Mekong River basin. Black square: Henicorhynchus siomensis. Star: H. lobotus. A total of 29 sites are serially 
numbered toward upstream within each of seven regions: N= Thai North, S= Songkhram, G= Gam, U= Mun, X= Xekong, T= Tonle Sap and M = 
IVlekong mainstream (MM). Location of the Khone Falls is indicated. 
doi:10.1371/journal.pone.0103722.g001 



sized cyprinids, collectively referred to as Siamese mud carp, are 
the most abundant and most economically important fish in the 
middle and lower Mekong basin [19], [20]. They are harvested in 
huge numbers, especially in and around Tonle Sap Lake in 
Cambodia and the Khone Falls area in southern Laos. These 
species account for 43% [1 1], [14] and over 50% [13] of the total 
catch in these areas, respectively, with an overall basin-wide catch 
being >12% for the two species combined [5]. 

The Siamese mud carp populations in the Mekong, or at least 
some of them, are known to perform long-distance migrations [4] , 
[7], [12]. However, what is known about their migrations is 



severely limited. In Thailand and Laos, the species undertake 
upstream migration in the early rainy season above the Khone 
Falls, whereas in Cambodia, this species migrate upstream at the 
onset of the dry season below the waterfalls [4]. The latter 
migrants possibly originate from the Tonle Sap Lake, first 
descending the Tonle Sap River and then ascending the Mekong 
mainstream toward and even past the Khone Falls [12]. Spawning 
takes place in tributaries or floodplains during the wet season, with 
a peak in May-June [21]. Eggs and larvae are carried to nursery 
habitats on the floodplain by the water current [4]. At the 
beginning of the dry season, juveniles move out of the floodplains 
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with the receding water and seek dry season refuge habitats such as 
deep pools in the Mekong mainstream. 

Multiple genetically discrete populations or stocks have been 
identified in both species of Siamese mud carp within the Mekong 
basin [22], [23]. However, ecological evidence that supports the 
genetic differentiation other than the single natural barrier of the 
Khone Falls has been lacking, thus hampering the development of 
effective fisheries management strategies. 

We used an otolith microchemistr^' technique to elucidate the 
natal origins and subsequent movement patterns of the Siamese 
mud carp in the Mekong River. Otoliths (ear bones) are paired 
calcified structures that continually grow by accreting calcium as 
well as other elements from ambient waters onto their surface [24] . 
Because of aceUular and metabolicaUy inert properties, the otoliths 
provide a chemical chronology over the entire life of a fish [25]. 
Although otolith microchc'mistry ^\-as originally dc\Tloped for 
marine or diadromous fishes that are exposed to higher elemental 
concentrations than freshwater fishes or undergo a substantial 
change in elemental concentrations during their life histories, the 
technique has increasingly been applied to obligate freshwater 
species. For example, to locate their natal origins [26-28] or 
maternal origins [29], to elucidate migration patterns [30-32], to 
discriminate wild from hatchery populations [33], and to detect 
exposure to pollutants from mine tailings [34]. Based on extensive 
field sampling throughout the lower Mekong basin, we compared 
the migration patterns of the Siamese mud carp between species as 
well as among geographical regions with particular emphasis on 
predicting the potential impacts of the proposed hydropower 
projects on these highly valuable fisheries resources. 

Materials and Methods 

Study area 

The study was conducted in the lower Mekong basin of 
Thailand, Laos and Cambodia (Fig. 1). The study area ranged in 
elevation from about 10 m a.s.l. at Kampong Chhnang, 
Cambodia (Tl) to 350 m a.s.l. at Chiang Saen, Thailand (N2). 
River discharge generally peaks in September, decreases toward 
the end of the dry season in April and starts to increase in May. 
Fish sampling was conducted during 2007-2010 at the total of 29 
sites in the mainstream and tributary systems in the following 
seven regions: the Mekong River in the Thai North (Nl-2), the 
Songkhram (Sl-4), Gam (Gl-4), Mun (Ul-8), and Xekong Rivers 
(XI -2), Tonle Sap Lake (Tl-5) and the Mekong mainstream in 
southern Laos and Cambodia (Ml -4, hereafter referred to as 
MM). 

Immediately below M3 in southern Laos near the Cambodian 
border are located the Khone Falls, the only natural waterfalls in 
the Mekong River that diverge into series of cataracts across the 
river width. Despite the large elevation gap of 21 m, one side 
channel called Hoo Sahong bypasses the waterfalls without any 
natural barriers along its approximately 7 km length, providing 
the single most important fish passage route year-round for the 
migratory fish species of the Mekong River [12], [35]. 

Most of the tributaries sampled were freely accessible from the 
Mekong River except for the Gam and Mun Rivers. In the Gam 
River, there are five irrigation dams located with approximately 
equal spacing along its 100-km long course before emptying into 
the Mekong River [36]. The most recent of the five dams was 
constructed in 2009 in the lowermost reach of the Gam, about 
2 km from the Mekong confluence. In the Mun River, a large 
hydroelectric dam, Pak Mun dam (17 m high; 136 MW), has 
blocked river water and possibly fish migrations at about 6 km 
from the Mekong confluence since 1994. The sluice gates of the 



dam have been kept open during the rainy season each year to 
allow for fish passage to mitigate adverse impacts on local fisheries 
[37], [38]. 

Several individuals of Henkorhynchus siamensis and H. lobalus 
were purchased at a fish landing site or the nearby local market, 
ensuring that the fish were captured in the vicinity and fresh (a few 
hours or less after landing). AU the fish were measured for fork 
length to the nearest millimeter and weighed to the nearest gram. 

Chemical analysis of otoliths 

A pair of otoliths (lapilli) of each fish was carefully extracted 
using ceramic forceps within a few hours of fish sampling. In 
Ostariophysi species such as cyprinids, lapilh are preferred otoliths 
to analyze because they are the largest of the three otoliths of a 
fish, have the single leading growth axis, and undergo conservative 
shape change throughout life [39] . The otoliths were immediately 
stored in a 2-ml polypropylene tube separately for each fish for 
transportation to the National Institute for Environmental Studies, 
Japan. In the laboratory, the otoliths were submerged into 30% 
hydrogen peroxide for 30 minutes and moved to a tube filled with 
ultrapure water (Milli-Q) for ultrasonic washing for 5 minutes to 
remove soft tissue around the otoliths. The otoliths were then triple 
rinsed with MiUi-Q, water and air dried in a clean booth equipped 
with HEPA filters. One of the pairs of lapilli was mounted whole 
on a glass slide with sulcus side down using thermoplastic glue 
(Crv'stalbond) for otolith surface assays. The other lapillus was 
embedded with high-transparency polyester resin (Technovit 
7200VLC, Heraeus Kulzer). After complete polymerization under 
fluorescent light for about 24 h, the embedded lapillus was sUced 
into a thin section in the frontal plane using a low-speed diamond 
wheel saw (SBT650, South Bay Technology). The sectioned 
sample was then mounted on a glass slide using Crystalbond for 
grinding with sandpaper and polishing with diamond slurry of 9 
and then 3 |tm to a thickness of 150-200 (im until the otolith core 
(i.e., primordium) was c-xposcd. The mounted sectioned sample 
was then sonicated, rinsed with Mill-Q_ to remove grinding 
particles and air dried. Both whole and sectioned otolith samples 
were kept in the clean booth until analysis. The sectioned samples 
were assayed only when greater than or equal to five otoliths were 
obtained at a given site-date combination. 

We used laser ablation inductively coupled plasma mass 
spectrometry (LA-ICP-MS; Agilent 7500cx ICP-MS, coupled 
with New Wave Research UP213 system equipped with a 
Nd:YAG 213-nm laser) to assay ^'^Na, ^^Mg, ''Ca, ''Mn, '"'"Zn, 

Sr, and ' Ba in the otolith samples. Otolith calcium (" Ca) was 
used as an internal reference to correct for any changes in 
materials ablated from sampling location to sampling location. A 
standard material produced by the National Institute of Standards 
and Technology (NIST 612) was used as an external reference to 
adjust for possible instrumental drift as explained below. Laser 
parameters were 80-100% laser power and 20 Hz repetition rate 
with a dwell time of 2 s per spot for both whole and sectioned 
samples. 

For whole otolith samples, a total of five spots of 90-nm 
diameter were ablated by directly applying the laser beam 
perpendicularly to the otolith surface to determine average 
elemental concentrations. For sectioned otolith samples, a series 
of 55-|.im diameter spots each separated by 60 |tm were ablated 
along the core-to-edge transect of the longest axis; on average 
24±3 (mean ± SD) spots were ablated depending on the otolith 
size. Data from ablation spots that overlapped, even partially, with 
surrounding embedding material near the otolith edge, were 
excluded from further analyses. Finally the otolith surface 
elemental data were combined with the transect data for each 
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otolith by adding the former to the last (i.e., edge) data point of the 
latter, producing a complete otolith profile of the fish'.s birth-to- 
death chemical history. Despite different ablation techniques 
adopted, the surface data were smoothly connected to the transect 
data without any noticeable gaps. 

The assay of a batch of five otoliths were preceded and followed 
by those of the NIST glass for calibration [40]. AH LA-ICP-MS 
analyses involved background data collection for 40 s, followed by 
otolith data collection for 60 s per spot. Helium gas stream was 
used to carry ablated materials from the laser system to the ICP- 
MS where it was mixed with Ar gas. Analytical accuracy of 
elements measured as the mean percent relative standard 
deviation (RSD) based on the NIST glass analyses were 8.3% 



(Na), 13.3% (Mg), 9.0% (Mn), 24.1% (Zn), 6.6% (Sr), and 7.3% 
(Ba). Limits of detection (LOD) were calculated for each element 
in each otolith sample as the blank value plus 3SD of the blank 
signals during the background data collection. Percentages of the 
samples that exceeded the LODs were 99.0/100 (Na), 100/100 
(Mg), 80.7/72.35 (Mn), 99.0/89.9 (Zn), 100/100 (Sr), and 100/ 
100 (Ba) for whole/sectioned samples. 

Water sampling and chemical analysis 

River water was sampled at a total of 74 sites throughout the 
lower Mekong basin by encompassing all regions but the Thai 
North where only otoliths were sampled (Fig. 2). Water was 
collected at the same or closest site to a fish sampling site when the 




Figure 2. Water sampling sites in six regions. S= Songkhram, G= Gam, U= Mun, X= Xekong, T= Tonle Sap and M= Mekong mainstream. 
doi:1 0.1 371/journal.pone.01 03722.g002 
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fish were bought far from the river. The sampling was conducted 
either from a boat or a bridge to avoid tlie direct influence of 
residential effluent alongside the river bank. Additional water 
samples were collected from up- and downstream of the original 
water sampling site when the confluences of large tributary systems 
existed in the vicinity to account for all potential water sources to 
which fish could have been exposed. The water samples were 
repeatedly collected in some regions to assess seasonal patterns in 
elemental signatures. Due to difficulties in accessing a river, 
however, samples could not be taken from some fish sampling sites 
in regions other than the Thai North. Conversely, the Mekong 
mainstream was extensively sampled for water from nearly as 
north as the Thai North to the Tonle Sap River confluence to the 
south. Immediately after the sampling, the water (14.5 ml) was 
filtered with 0.20 |im syringe filter into an acid-washed 15 ml 
polystyrene test tube to remove microorganisms and particles, 
acidified with 0.5 ml of aristar grade nitric acid to 1% v/v, and 
refrigerated until analysis. 



The water samples were analyzed for the concentrations of the 
same set of elements as otoliths using both an Inductively Coupled 
Plasma Atomic Emission Spectrometer (ICP-AES) (ICAP-750, 
Nippon JerreU-Ash Co. Ltd.) primarily for major and minor 
elements and an ICP-MS (Agilent 7500cx) for trace elements. 
Some elements were assayed by both instruments, but when ICP- 
MS readings exceeded 30 ppb, corresponding values from ICP- 
AES were adopted. The detection limits (ppb) of ICP-AES were 
100 (Na), 20 (Mg), and 30 (Ca), and those of ICP-MS were 0.001 
(Mn), 0.004 (Zn), 0.001 (Sr), and 0.003 (Ba) with all sample 
readings being above the detection limits. Elemental concentra- 
tions both in water and otoliths were converted to molar ratios of 
analyte/Ca in all statistical analyses. 

Data analyses 

Molar ratios of Mg, Mn, and Zn to Ca for both water and 
otolith samples were logarithmically transformed to achieve 
normality before statistical analyses. Na was excluded from all 
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Figure 3. Sr-Ba elemental signatures of river water samples collected from six regions during 2007-2010. 
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the analyses except water-otolith linear regression analysis because 
of its significant multicoUinearity with most of the other elements. 

To eliminate random noise from the otolith profile data of 
elements from the core to surface, the Gaussian-type 5-point 
running weighted average was applied to smooth the original data 
with a weight being equal to: 

G(x)=exp(-^)^^|:^exp(-0 

where an integer x varied from —2 to 2. This resulted in the losses 
of 2 consecutive original data points (or 120 |lm) each from the 
core and edge portions of the profile, although the lost data were 
partially reflected in the first and second data points of the core 
portion as well as the last and second last points of the edge portion 
in the resultant smoothed profile. 

Water sampling sites located within a 1 0-km radius of each fish 
sampling site were identified using the geographical information 
system (ArcGIS, ESRI) to calculate average water elemental 
concentrations for the fish site. Linear regression analyses were 
performed to examine relationships between water and otolith 
surface elements. The effects offish size and sampling season on 
the water-otolith relationships were tested by incorporating fish 
size to the linear regression as a candidate explanatory variable 
and by constructing nested regression models with separate slopes 
and intercepts between wet (June-November) and dry (December- 
May) seasons, respectively. A partition coefficient of elements was 
estimated as a slope of a regression line through the origin. 

Regional and temporal differences in water and otolith surface 
elemental signatures were tested with the two-way analysis of 
variance (ANOVA). Multivariate analysis of variance (MANOVA) 
was also performed with PiUai's trace statistics to compare otolith 
elemental signatures among regions and between species. Multiple 
comparisons among regions were made for selected elements using 
the Tukey's HSD test separately for otolith core and surface. 
Linear discriminant analyses (LDA) were applied separately to the 
water and otolith data to examine classification performance of the 
sampling regions. For the otolith LDA, only H. siamensis samples 
from aU regions but the Xekong River, were used due to the 
hmited samples available from this region. Elemental molar ratios 
were standardized to have the mean 0 and standard deviation 1 to 
allow the comparison of relative contributions by each element in 
the LDA. Jackknife cross-validation procedure determined classi- 
fication accuracy of the regions. All the statistical analyses were 
conducted with a significance level of a = 0.05 using S-PLUS 
Version 8.1 (TIBCO Software Inc.). 

Ethics statement 

Field sampling was carried out with the relevant permissions 
from Department of Fisheries, Thailand, Liland Fisheries 
Research and Development Listitute, Cambodia, and Living 
Aquatic Resources Research Centre, Laos. This study did not 
involve any endangered or protected species. The species involved 
were originally collected by fishermen for human consumption 
using the approved fishing gear in respective countries. Therefore, 
no approval was obtained from any animal ethics committees. 

Results 

Water elemental signatures 

A total of 158 water samples were collected throughout the 
lower Mekong basin (Table SI). Water elemental signatures 
dififered among regions in each of five elements (i.e., Mg, Mn, Zn, 
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Sr and Ba; p<0.05) and among years in Mg, Sr and Ba (p<0.05). 
To a large extent, the six regions were separated with two 
elements, i.e., Sr and Ba (Fig. 3). The LDA of the two-element 
model classified 68.4% of the data c:orrecdy to the regions of 
sampling, whereas a full model based on the five elements 
classified 86.7% of the data (Table 1). The correct classification 
rates were particularly high for the Mekong mainstream (100%) 
despite its extensive sampling area and the Mun River (98%), for 
which the percentages were not different between the models. 
Although the Songkhram River appears to have shifted its water 
Sr-Ba signatures over the years, the other regions were relatively 
consistent in the chemical signatures. 

Water-otolith relationships 

A total of 1 73 Henicorhynchus siamensis and 24 H. lohalus were 
collected from the seven regions of the Mekong basin (Table S2). 
Fish fork length varied both among regions and sampling months 
with the overall mean ± SD of 131 ±22 mm for H. siamensis and 
130±15 mm for H. lohatus. Water samples collected within a 10- 
km radius of fish sampling sites were available for 1 22 H. siamensis 
and all the 24 H. lohalus otolith samples. Linear regression lines 
between water and H. siamensis otolith surface elements were 
significant and positive for Zn, Sr and Ba, while for H. lobatus the 
only significant and positive regression was found for Sr (Fig. 4). 
Two regression lines for Sr were not different between the species 
both in intercepts and slopes (p>0.05). Partition coefficients were 
0.043 for Zn (p = 0.001), 0.409 for Sr (p<0.001) and 0.055 for Ba 
(p<0.001) for H. siamensis, and 0.496 for Sr (p<0.001) for H. 



lohatus. Fish size had no significant effect on any of the water- 
otolith elemental relationships (p>0.05). However, the sampling 
season had a significant effect on Zn and Ba, such that only the 

otoliths and river water collected during the wet season (i.e., June- 
November) had significantiy positive correlations (p = 0.001 and < 
0.001, respectively). 

Otolith elemental profiles 

Of the five elements in H. siamensis otolith surfaces, Mg, Mn, Sr 
and Ba were difl'erent among regions (Mg: 160 = 5.183, p< 
0.001; Mn: F^^ 16,, = 5.864, p<0.001, Sr: 'igo = 27.300, p< 
0.001; Ba: F^^ 'i6o = 4.353, p = 0.001), and Mn, Sr and Ba were 
difierent among years as well (Mn: i6o = 2.942, p = 0.035, Sr: 
^3, 160 = 7.871, p<0.001, Ba: F3, i6o = 6.352, p<0.001). The 
classification accuracy of the two-element (Sr and Ba) LDA model 
was low (39. 1 %), and that of the four-element model was still not 
very high (50.3%, Table 2). Note, however, that 17 out of 21 
otoliths (81%) from the Gam River were correctiy classified. 

Otolith Sr and Ba for H. siamensis varied largely throughout 
life, but the shape of the elemental profiles resembled one another 
among indi\'iduals in a gi\'en site and even a region (Fig. 5). 
Between regions, however, otolith surface elemental signatures 
were different (Pillai's trace =0.814, f 10, 202= 13-872, p<0.001), 
and core signatures were also different (Pillai's trace =0.765, 
F\a, 202 =12.520, p<().001). Otolith surface Sr signatures were 
different in all pairs of regions involving the Thai North and Gam 
River (Table 3). fn particular, the Gam River exhibited the lowest 
level of otolith surface Sr among all regions. Otolith surface Ba was 
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Figure 4. Otolith surface elemental signatures plotted against river water signatures. Open circle: H. siamensis. Closed circle: H. lobatus. 
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similarly different in many pairs involving these two regions. The 
otolith core Sr level was significantly lower again in the Gam River 
than the other regions (Table 4), and to a lesser extent the Ba level 
was also low in this region. The otolith core Ba was highest in the 
Thai North among all regions. 

Fish from S 1 exhibited little variations in both elements possibly 
due to their small body size (<100 mm). Despite large body sizes 
(& 130 mm), however, the Gam River fish showed severely limited 
variations in the elements throughout their lives, the profiles of 
which were confined in small areas at very low Sr and Ba levels. 
Profiles for the Mun River, Tonle Sap and MM fish were similar 
to each other and generally increased both in Sr and Ba with fish 
growth. 

Among three sampling sites along MM (M2 below Khone Falls, 
and M3 and M4 above the falls), the otolith surface signatures of//. 
lohatus were diflFerent (PiUai's trace =0.811, 28 = 4.779, 
p = 0.005), but their core signatures were not (Pillai's trace = 0. 190, 
28 = 0.7 3 3, p>0.05. Fig. 6). Between the two species both 
collected at M2 in June 2010 (Figs. 5 & 6), otolith surface signatures 
were not diflFerent (Pillai's trace = 0.459, f 2, 9 = 3.824, p>0.05), but 
core signatures were diflFerent (Pillai's trace =0.536, ^2, 9 = 5.196, 
p = 0.032). Most profiles in bodi species commonly exhibited a slight 
decrease (ca. 0.2 mmol/mol) in Sr/Ca immediately outside of 
primordia at a very early life stage of fish, which was followed by 
increases in Sr/Ca and Br/ Ca at later life stages. 

Discussion 

Uptake of elements by Siamese mud carp otoliths 

The Siamese mud carp undertook complex migrations as 
revealed by Sr-Ba profiles along the otolith growth axis. There 
were similarities in the profiles among individuals within regions 
and dissimilarities between regions, suggesting that fish from a 
given region hatched out in a similar environment and migrated in 
schools before being harvested by fishermen, but that their natal 
origins and migration pathways were diflFerent across most regions 
throughout the Mekong basin. 

Similarly to many previous studies, strontium was the most 
reliable marker element followed by barium to elucidate Siamese 
mud carp migrations in the Mekong River [26], [30], [41], [42]. 
Strontium partition coefficients for H. siamensis (0.409) and H. 
lohatus (0.496) were comparable to those observed for freshwater 
fishes in temperate and arctic climates such as cutthroat trout 
Oncorhynchus clarki lewisi (0.40) [26] and Arctic grayling 
Thymallus arcticus (0.369) [30]. Barium partition coeflFicient for 
H. siamensis (0.055) was again similar to those for cutthroat trout 
(0.04) and Arctic grayling (0.051). Magnesium incorporation to 
otoliths was not significantly correlated with water Mg/Ca. 
However, the negative, although insignificant, relationship be- 
tween water and otolith Mg we found in H. siamensis was also 
observed in the otoliths (and scales) of cutthroat trout [26], for 
which an underlying mechanism has yet to be understood. We 
found no correlation between water and otolith Mn, which was 
similar to the findings in previous studies [41-44]. Although a 
mechanism for Mn incorporation into otoliths is largely unknown, 
physiological processes such as dietary uptake have been 
implicated [42], [43]. 

A temperature effect on otolith microchemistry is generally less 
important than ambient elemental concentration and depends on 
the ranges of both temperatures and elemental concentrations 
[44-47] . Mugiya and Tanaka [48] observed a positive influence of 
water temperature on otolith Sr uptake in goldfish (cyprinid), but it 
was on the order of one fifth of Sr/ Ca in ambient water. Miller 
[47] observed a positive temperature effect on Sr (and Mg) uptake 
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133±4 
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Figure 5. Sr-Ba elemental profiles along H. siamensis otolith transect. Triangle: otolith core. Circle: otolith surface. From top to bottom, 
otolith samples from Thai North (2), Songkhram (3), Gam (2), Mun (2), Tonle Sap (3), and Mekong mainstream (2) were analyzed with the number of 
sites in a parenthesis. Within a given site, different colors denote different individuals. 
doi:10.1371/journal.pone.0103722.g005 

Table 3. Pairwise differences {mmol/mol) between regions in H. siamensis otolith surface elemental signatures for Sr/Ca (above 
diagonal) and Ba/Ca (below diagonal). 







Thai North 


Songkhram 


Gam 


Mun 


Tonle Sap 


MM 


Thai North 




-0.499 


0.448 


-0.716 


-0.557 


-0.395 


Songkhram 


-0.033 




0.948 


-0.217 


-0.058 


0.104 


Gam 


-0.055 


-0.022 




-1.16 


-1.01 


-0.844 


Mun 


-0.021 


0.012 


0.034 




0.159 


0.321 


Tonle Sap 


-0.038 


-0.005 


0.017 


-0.017 




0.162 


MM 


-0.029 


0.004 


0.026 


-0.008 


0.009 





Significant differences are in bold. Values are average elemental molar ratios for regions in the row minus those in the column headings. 
doi:l 0.1 371 /journal.pone.Ol 03722.t003 
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Table 4. Pairwise differences {mmol/mol) between regions in H. siamensis otolith core elemental signatures for Sr/Ca (above 
diagonal) and Ba/Ca (below diagonal). 







Thai North 


Songkhram 


Gam 


Mun 


Tonle Sap 


MM 


Thai North 




0.141 


0.757 


0.063 


-0.071 


0.166 


Songkhram 


-0.009 




0.616 


-0.079 


-0.212 


0.024 


Gam 


0.030 


-0.022 




0.694 


0.828 


0.592 


Mun 


0.025 


-0.016 


0.005 




-0.134 


0.103 


Tonle Sap 


-0.021 


-0.012 


0.009 


0.004 




0.236 


MM 


0.038 


0.029 


-0.008 


-0.013 


-0.017 





Significant differences are in bold. Values are average elemental molar ratios for regions in the row minus those in the column headings. 
doi:l 0.1 371/journal.pone.Ol 03722.t004 



in juvenile Chinook salmon Oncorhynchus Uhawytscha otoliths, 
but the effect was observed in the Sr-rich saline environment and 
not in fresh water. 

Migrations by region 

Despite some limitations of and uncertainties about the 
technique, we could detect region-specific as well as common 
basin-wide patterns in the migrations of Siamese mud carp. In the 
Thai North, H. siamensis exhibited greater variability in Ba rather 
than Sr throughout their life history, resulting in a peculiar shape 
in the elemental profiles. This may be explained by a diflFerence in 
geological structure between upper and middle-to-lower Mekong 
basins. The upper basin including the Thai North is characterized 
by multiple tectonic plates that are intricately separated by several 
major Tertiary-Quaternary sedimentary basins, whereas the 
middle Mekong basin, especially the vast area called the Khorat 
Plateau, encompasses three regions, i.e., Songkliram, Gam and 
Mun river basins in northeastern Thailand, and is covered 
uniformly by Cretaceous sedimentary rocks [7], [49]. Therefore, 
the larger otolith Ba variability does not necessarily translate into 
longer fish migrations but reflects a greater spatial complexity in 
water elements encountered by the fish. Due to the unique otolith 
elemental profiles, H. siamensis from the Thai North likely did not 
share migration routes with individuals from the other regions, 
corroborating a finding by Adamson et al. [23] who identified a 
genetically distinct population of the species in this region. 

Otolith core as well as surface elemental signatures were not 
different between the Songkhram and Mun Rivers. Adamson et al. 
[23] considered H. siamensis from these regions as a single 
genetically homogeneous stock. Although both the mean otolith 
chemical signatures and the genetic analysis could not differentiate 
between the two regions, the shape of the otolith elemental profiles 
was considerably different with much more diverse core signatures 



for the Songkhram fish in the dry season than the Mun fish, 
suggesting that multiple H. siamensis stocks utilize the Songkhram 
River in this season. Located between the Songkhram and Mun 
Rivers, the Gam River had a H. siamensis population with 
extremely low levels of, and a limited variability, in otolith Sr and 
Ba during their life history. This will be discussed later in the 
context of a possible damming influence on this population. 
Elemental profiles largely overlapped between Mun, Tonle Sap 
and MM with no significant differences in Sr or Ba at both otolith 
core and surface. However, because the genetic analysis identified 
a single stock of//, siamensis below the Khone Falls [23], H. 
siamensis in the Mun River, which connects with MM above the 
Khone Falls, likely have a different natal origin from fish in Tonle 
Sap or MM. During the peak flood season, the Mekong River 
water causes the Tonle Sap River to reverse its flow [50], likely 
mixing the water and homogenizing elemental signatures 
throughout the vast geographical extent. This unique hydrological 
event may have contributed to the similar otolith elemental 
profiles between the Tonle Sap and MM regions. However, it does 
not preclude the possibility of the species collected in MM 
originating from, or spending part of their life in, the Tonle Sap 
Lake [11], [12]. 

H. lobatus collected at M2 (ca. 70 km below the Khone Falls), 
M3 (right above the falls), and M4 (ca. 130 km above the falls) 
exhibited fairly synchronized indistinguishable Sr-Ba profiles with 
the core elemental signatures being statistically not different 
among the sites. This finding supports the previous notion that, 
unlike H. siamensis, H. lobatus can traverse the Khone Falls [12] 
and that no genetic differentiation exists among individuals 
collected below and above the waterfalls [22]. Although the 
similarity in the otolith profiles and particularly core elemental 
signatures could be explained by the passive drift of H. lobatus 
from above the Khone Falls downstream to M2 below the 
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Figure 6. Sr-Ba elemental profiles along H. lobatus otoWVn transect. Triangle: otolith core. Circle: otolith surface. All samples were collected in 
June 2010 from 3 sites in the Mekong mainstream. Within a given site, different colors denote different individuals. 
doi:1 0.1 371/journal.pone.01 03722.g006 
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waterfalls, Hurwood et al. [22] considers it unlikely because the 
haplotypic diversity of the species was higher above the Khone 
Falls than below the falls, indicating the active upstream dispersal 
of the species over the waterfalls. In summary, the Khone Falls are 
probably not a migration barrier at least to H. lobatus as local 
ecological knowledge [4], catch statistics [12], genetic analysis [22] 
and otolith microchemistry have all indicated. 

Otolith elemental profiles of the two species commonly 
exhibited a sharp drop in Sr level just outside of a primordium 
at a vcr)- early life stage. It spanned approximately 20()-3()() |j,m of 
an otolith distance from the core or er[ui\'alently l.^-2()'/() of the 
average fish life as an average otolith core-edge length was 
1,440 ^lm. Although it only takes 12-16 hrs for Siamese mud carp 
to hatch out after spawning [51], it takes seven days for their fin 
rays to start developing and 35^0 days to fully develop [52]. 
During this early developmental stage, they would be flushed 
downstream by water currents for long distances before becoming 
fuUy capable of swimming against river flows. However, relating 
the initial drops in otohth Sr to passive movement of fish with 
water currents may not be a straightforward explanation because 
water Sr/Ca does not necessarily decrease in a longitudinal, 
downstream direction in the Mekong River and its tributaries [53] . 
Nothing is known about finer-scale lateral variations in Sr/Ca 
across a river channel and floodplain where the species are 
believed to rear, so that the passive drift of fish eggs and larvae 
cannot be entirely ruled out as the possible explanation of the 
phenomenon. Alternatively, it could have resulted from reduced 
Sr uptake by the otoliths of Siamese mud carp larvae at the onset 
of the rainy season when water temperatures start to drop. An 
annual temperature variation is roughly eight degrees in Celsius in 
the Mekong tributaries (M. Fukushima, unpublished data), which 
is substantial for tropical ecosystems. Despite relatively minor 
influence of temperature on otolith elemented uptake as described 
above, a rapid decrease in water temperature during the early 
developmental stage may have played a role. 

We found that Zn and Ba had insignificant correlations between 
water and otolith surface during the dry season. It may indicate 
more active, faster migration by the species in this season, as fish 
would spend less time at any given location for elements unique to 
that location to be incorporated into otoliths. Otolith trace 
elements such as Sr and Ba change rapidly in response to changes 
in water chemistry, but it t\7)ically takes about two weeks for the 
otolith concentrations to stabilize and fuUy reflect ambient 
concentrations [48], [54], [55]. 

Effects of damming on Siamese mud carp populations 

The migration of H. siamensis was severely suppressed in the 
Gam River. OtoUth profiles for this region not only occupied a 
small confined area in the Sr-Ba biplots but also overlapped 
minimally with profiles for the other regions including MM. Fish 
sampling sites in the Gam River were located in-between irrigation 
dams that must have severely impeded fish movement within the 
tributary as well as movement to MM. Although a fishladder has 
been installed at each dam site in this river, its performance has 
received a low rating because of its limited efficacy favoring only 
small-sized short to medium-distance fish migrants [36]. H. 
siamensis are known to be able to establish populations in a lentic 
system such as a reservoir [56] and even river reaches between 
barriers [4] as they probably do in the Gam River judging from 
the large body size of individuals collected. Because of this 
flexibility, the species has been among the candidates for fish 
stocking to enhance fisheries production in reservoirs in Thailand 
and Laos [57]. However, their gonadal as well as somatic growth is 
significandy slower than those observed in wild populations in lotic 



systems [56], [58], and therefore they do not flourish in those man- 
made water bodies [59], [60]. 

In contrast, H. siamensis migration appeared not to be 
suppressed in the Mun River despite a possible impediment by 
the Pak Mun dam near the confluence of the Mekong River. 
Water Sr-Ba signatures in the Mun River and MM were relatively 
constant throughout a year and across years within the regions, 
and they were fairly distinct from each other. Nevertheless, the 
otolith Sr-Ba profiles of the two regions considerably varied and 
overlapped, suggesting that the Mun fish had migrated as 
extensively as did MM fish through the sluice gates of the Pak 
Mun hydroelectric dam, presumably from MM [58]. This study 
reaffirmed the importance of opening the sluice gates during the 
flood season to conserve fish species diversity and local fisheries 
[37], [38]. 

H. lobatus win suffer from the planned Don Sahong 
hydropower project in the middle of Hoo Sahong channel, the 
single most important fish migration corridor bypassing the Khone 
Falls. Otolith Sr-Ba profiles indicated little difference in the natal 
origin and subsequent migration pathways of the fish collected 
below and ab()\'e Khone Falls. Along with H. lobatus, other 
commercially important fish species would be equally affected by 
this dam. Pangasius krempfi, highly valued large-bodied fish 
species mostly harvested around the Khone Falls area during their 
spawning migration, were recendy confirmed by otolith micro- 
chemistr\' to be anadromous, down-migrating to tlu; s('a for 
feeding and up-migrating the Mekong River over the Khone Falls 
for reproduction [61], Artificial fish passes are unlikely to be an 
effective mitigation measure in the Mekong River because of the 
size offish migrations which can reach as much as 30 tonnes per 
hour [62] and because of the significant number of migratory fish 
species that differ in body size and shape and behavior [5]. 

The impact of this and many other hydropower projects 
proposed in the Mekong basin on Siamese mud carp populations 
can be devastating in the light of the findings of this study. First, 
the species have discrete and rigid migration pathways that are 
consistent among individuals within a region but different among 
regions, making it unlikely that fish seek alternative migration 
routes once a dam blocks their original migration route. Second, 
their migration range has already been severely suppressed in a 
dammed tributary of the Mekong, i.e., the Gam River, despite the 
existence of fish passes. Third, there is a population of H. lobalus 
distributed along the >200-km stretch of the Alekong mainstream 
that are migrating over the Khone Falls, most likely through the 
Hoo Sahong channel where the Don Sahong dam is planned. 
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